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Abstract 

For a substance diffusing on a curved surface, we obtain an explicit relation valid for very 
small values of the time, between the local concentration, the diffusion coefficient, the intrinsic 
spatial curvature and the time. We recover the known solution of Pick's law of diffusion in 
the flat space limit. In the biological context, this result would be useful in understanding the 
variations in the diffusion rates of integral proteins and other molecules on membranes. 
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1. Introduction 

Transport of enzymes, charged ions and metabolic substances within biological cells and tissues and 
across cell membranes is one of the major processes which sustains and guides life. Indeed, extra- 
cellular and intracellular transport of substances can well be considered to be the most important 
and pervasive among all the life-supporting biological activities. 

Molecular transport across cell membranes by passive diffusion or in accordance with Pick's law 
is a well-studied area. However in the available literature on the subject, no mention has been made 
of how the local curvature of the cell plays a part, if at all, in this process. Molecules released at a 
specific location on the cell surface or on the nuclear membrane diffuse along the curved membrane 
surface to another location. 

It is known that thermal agitation permits lateral diffusion of phospholipid and glycolipid molecules 
within a leaflet of planar phospholipid bilayers of biological membranes. A lipid molecule can diffuse 
several micrometres per second at a temperature of 37 degrees C. It has also been established 
experimentally that many important proteins freely float within the plane of the membrane. 
Measurements have shown that the rates of diffusion of proteins in biomembranes are considerably 
lower than those seen in artificial membranes [1,2]. The physical structure and the dynamical changes 
occurring on a membrane surface would well be expected to play an important role in determining 
the lateral mobility of molecules on its surface. 

The metabolism and synthesis of fatty acids and phospholipids occur in the smooth endoplasmic 
reticulum, and the rough endoplasmic reticulum is a site of protein synthesis. It is well known 
that in many cells these extensively curved and folded membrane vesicles are continuous with the 
nuclear and cell membranes. In the c}rtosol also, these folds distort the homogeneity in the spatial 
distribution of the cytosolic fluid. 

Transport of a substance by diflFusion should therefore be described by a corrected form of Pick's 
law, modifled to take into account the local curvature of the surface through which it moves. 

In this paper we discuss how to take care of curvature effects and also give for transient phenom- 
ena, an explicit expression relating the concentration of the diffusing substance, the intrinsic spatial 
curvature experienced by it, the diffusion coefficient and the time. 
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2. Diffusion on curved surfaces 

Consider diffusion of a substance described by its concentration C{x, t) from a spatial point x where 
it has been released on the cell, to another point x' . For a particular time slice, the line element ds 
between each pair of neighbouring points on the spatial surface is given by: 



ds^ = ^ gij{x)dx'dx^ (1) 

where dx'^ denote the coordinate differences between neighbouring points, n is the spatial dimen- 
sion and Qij denotes the metric. We choose to work with a Riemannian signature for the metric. 
The usual form of Pick's law relates the current density or the flux of material per unit area, j{x, t) 
to its concentration gradient in flat space : 

ji{x,t) = -DdiC{x,t) (2) 

where D denotes the diffusion coefScient, 9, denotes the gradient operator, and C{x, t) is the field 
variable denoting the concentration. It is assumed here that the diffusion coefficient is independent 
of the concentration of the diffusing substance. 

In curved space, while formulating the problem, one must incorporate the effects of the intrinsic 
spatial curvature of the surface on which the substance is diffusing. 

We make the simplifying assumption that in the infinitesimal neighbourhood of any point, the 
diffusion properties are the same in all directions and that D does not depend upon the position 
and the concentration of the diffusing material. 

Transport of the substance by diffusion into and out of the invariant volume element 
surrounding the point x is given by the conservation equation : 

^ = -V./(x,*) (3) 
where Vi denotes the covariant derivative and includes the Christoffel connection T'^u : 

and we have considered a parametric dependence of C on the time t. 

Performing a covariant differentiation of (2) with respect to x and substituting (3) in it we then get 
the correct form of Pick's second law of diffusion : 

= -DOCix,t) (4) 
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where we have used □ to denote the n-dimensional Laplace-Beltrami operator. For flat 3-diniensional 
space, □ reduces to the usual 3- dimensional Laplacian. 

It has been shown in [3] that at least in mitochondrial inner membranes, the diffusion coeflacient D 
of intramembrane particles shows an inverse correlation with their concentration, implying that the 
proper form of Pick's law reflecting the concentration dependence of D should be studied, rather 
than equation (4). In our work however, we consider only the simplest form of Pick's law with a 
concentration-independent diffusion coeflacient, in order to see how far just the spatial curvature 
effects could modify the known result. 

It becomes particularly interesting to learn about the configuration of the released substance during 
the initial infinitesimal time intervals to see how the intrinsic curvature of the cell would influence 
diffusive transport on the membrane surface and hence its configuration at later times. 
In order to solve (4), we rescale the time parameter by: t ^ Dt so that (4) now reads: 

^£^ = DCix,t)-r,C{x,t) (5) 

and the parameter t now has the dimensions of length squared. We have introduced a drag term 
r/C with ?7 > which can be thought to account for negative concentration changes due to possible 
frictional effects on the motion of the molecules. We have introduced it here just for the sake of 
mathematical convenience and at the end of the calculations it can be set to zero. 
In the actual physical situation, of course, the drag term is very much present and gets contributions 
from the drag arising from the pericellular matrix viscosity, from steric effects, and from transient 
binding to relatively immobile structures [1]. Also in the actual situation, the t] term is not constant 
and has a coupling with the concentration gradient. We have however restricted ourselves to jj = 
for the sake of simplicity in this paper. 

We assume that the molecules diffuse freely on the surface without interacting or binding with any 
other molecules. We write equation (5) in a point-separated form as : 

dC{x,x',t) , 

= [U^ -'q)C{x,x ,t) (6) 

where the biscalar C{x,x',t) is subject to the condition : 



lim C{x,x',t) =C{x,t) (7) 

x'— >0 



and to the physical boundary condition: 



limC{x,x',t) =S{x,x') (8) 
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This enables us to get a well-defined explicit solution for C{x,t) which is valid for small values 
of t, in terms of the spatial curvature. The solution to (6) is well known [4-6]: 



C{x,x',t) = ^e-''*e-^^Ai/2(a;,a;')0(a;,x',i) (9) 

whore the biscalar (j{x, x') equals half the square of the geodesic distance between x and x' and 
lS.{x,x') is the VanVleck-Morette determinant : 

A(x,x') = -{g{x))-''^Aei[-d,dra{x,x')]{g{x'))-''^ (10) 
This is a biscalar quantity which reduces to unity in flat space. 

In curved space, one can expand A(a;, x') in a series expansion in powers of the curvature by working 
in Riemann normal coordinates y which define a locally inertial system in the neighbourhood of the 
point x' . In these coordinates [7], with origin at x' , IS.{x,x') = {g{x))~^/^ , so that 

Ay\x,x') = {g{x))-'/' = 1 + ^Ra^y^y^ + 0{y') (11) 

where x is regarded as a function of the Riemann normal coordinates y, such that X ^ 3C clS 
y — > 0. In the coincidence limit, and for our purposes, it is only the first term on the right hand 
side of (11) which is relevant for the calculations. 

The function Q{x,x',t) has the following series expansion in the coincidence limit x' ^ x : 

oo 

lim n{x,x',t) = y^t''Ek{x) (12) 

x' — *x ^ — ^ 

valid in the limit i — + where Ej. {x) are known coefficients known in the literature as Gilkey coefii- 
cients [4-7] : 



Eq = I 

R 

Es = ^ [iSa^R + nR.^.K'" - 2R^,,pRi"''P - 4i?^,.,pi?''^'^ + QR^^p^.rR""'"'''^ + 28RaR 

-8R^,^R^^'' + 2ARf,,R''P'% + 12i?^,^,ni?''^^^ + - ^RR^,R>^- + 1^RR^^^^R^^-P- 

- R P'^P -\-—R R f}W^^ -—R Rl^ nvpuT I D nhtvaP npa 

+fR^.paR>''''^R''a%] - ^o'v + ^v^v + ^{^v)v + ^V;^v'' - Iv' - ^R^v 
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Here, / denotes the identity matrix, R stands for the Ricci scalar and the semicolon denotes a covari- 
ant differentiation. Although the fourth Gilkey coefficient has also been calculated in the literature, 
we have displayed above only terms upto third order in the Riemann curvature. 
Rescaling now t back to Dt we obtain the solution we seek for diffusion of molecules in the presence 
of a drag term in n spatial dimensions for transient times : 



For the standard diffusion equation without the r] term, and in flat space, the only Gilkey coef- 
ficient which contributes is Eq and in this case we recover the known result: 



In a recent paper [8], Gompper and Goos have suggested that the diffusion of amphiphilic 
molecules within a monolayer at the oil-water interface of the microcmulsion phase in an oil-water- 
amphiphile mixture can be used to measure the average Gaussian curvature of the monolayer. They 
have considered surfaces of constant curvature. The result (14) we have discussed here for the 
concentration, is valid for surfaces of varying curvature also. 

In fact, from (14) it is an easy matter to get a general solution to the diffusion equation for 
n-dimensional spaces with arbitrary constant curvature K for which the value of the Riemann cur- 
vature depends neither on the coordinate x nor on the planar direction at x. For such spaces, the 
Riemann curvature is given in terms of their metric gij by 

Rtjki = K{g^k9ji - gugjk) (for n > 3) , (16) 
whence the Gilkey coefficients turn out to be : 



^0 = 1 

^ n(n-l)^ 
6 

^ ^ n(n- l)(3n-H) ^^2 
^ 360 

E3 = '^^^^^{7{n-lf{5n-l) + 61{n-lf + 68n + 28}K^ (17) 

We then find that the expression (14) for the concentration of the diffusing substance has the fol- 
lowing dependance on the Gaussian curvature K : 
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C(x,t) = j:re~^(l + -— '-KDt+— '^^ ' (KDty 

{AiiDtT'^ V 6 360 ^ ^ 

+ ^:^^^^[7(n-lf(5n-l) + 61(n-l)' + 68n + 28](i^£>t)' + ...) (18) 

K > Q corresponds to the spherical surfaces while surfaces with K < Q correspond to hyperboloid 
ones — the K = {) are flat Enclidean surfaces. 

It is shown in [8] that the structure of a microernulsion can be quantified in terms of a quantity 
which depends upon the Euler characteristic xe of the surface within which the amphiphile molecules 
diffuse. Xe is obtained from the Gaussian curvature using the Gauss-Bonnet theorem : 



/ 



dSK = 2-kxe (19) 



where the integral is over a closed siirface S. It should be therefore possible in the case of surfaces 
of approximately constant area, to express the result (18) in terms of the topological invariants 
characterizing them, after appropriately scaling them. This exercise is however beyond the scope of 
this report. Because of their enormous complexity, biological cells and membranes do not in general 
have isotropic and homogeneous composition, and the membrane surfaces arc more often than not, 
of varying curvature, and in these situations, one needs to use (13) and (14) rather than (18). 

For the specific case of diffusion in two dimensions such as on membranes, the coefficients in 
(13) simplify considerably because in these dimensions both the Riemann tensor Rijki and the Ricci 
curvature scalar R have only one component and both the Riemann tensor and the Ricci tensor i?^ 
can be expressed in terms of the curvature scalar R : 

Rijkl = ^R{9ik9ji - 9ii9jk) (20) 

and 

= 2^9ij (21) 

It must be borne in mind that in our treatment, we have regarded time as a parameter and the 
indices i, j, k, I, fi, v, a, f3 etc. label spatial dimensions only, for we are working on a particular time 
slice at each instant of time. 

We consider the simplest example of diffusion of a substance on the surface of a sphere of constant 
radius r. For a 2-sphere, the Ricci curvature scalar is: 
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In fact, in this case, the Gaussian curvature K = Substituting this value of R in (13), (20) and 
(21), or by simply using (17), the Gilkey coefficients reduce for this example to: 



Eo = 1 

1 



El 
E2 



3r2 
1 



15r4 

= 3lfc^ (22) 

giving the following result for the concentration of the diffusing substance of an initial unit amount, 
at a point distant x from the point of its release on the surface of the sphere, at a time t: 

''^'^'^ = j4^f~^i'+3^ + r5^^^ +315(7^) +•••) (23) 
The result obtained in [8] for the mean square displacement of a particle diffusing on a sphere is 

essentially equivalent to the leading and next to leading order terms in (23). 

In the series expansion in (14), (18) and (23), valid for very small t values, it is assumed that 

the curvature terms are small in comparison with the flat space result. Care must be taken before 

applying the actual values of t,D and r to these expressions to ensure that this assumption is 

satisfied. 

For a substance having a D value of 10~^cm2s~^ released on the surface of a spherical cell of 
radius 1/um , diffusing through a distance of 0.5/im in time 1ms, one obtains a calculated value of 
4259.4751 x 10^ per cm^ for its concentration, using the usual expression (15) for flat space diffusion, 
while the improved solution (23) gives an additional correction of 0.034 per cm^ to this - a difference 
of 3.4% from the flat space result, and a deviation of 0.33% from the flat space results for a time 
duration of 10~^s, while for a D value of 10~^cm2s~^, the deviations from the flat space result for 
time durations of 1ms and 0.1ms are 0.33% and 0.033% respectively. 

The experimentally measured values of the concentrations, of course, correspond to the corrected 
values and the curvature-corrected Pick's law (14), since the diffusing molecules have already tra- 
versed over the curved surface of the cell. However, it must be borne in mind that the diffusion 
times measured and calculated, in fact, calibrate distances different from the flat space distances, 
when one makes a comparison between diffusion rates on different cells and on membranes whose 
curvatures differ from point to point and from one another. What we intend to point out here is 
that one must remember that flat space methods must not be applied when one is talking about 
biological membranes and surfaces which are very curved, or even for surfaces with varying curva- 
ture for which one must apply the coefficients in (13). We have not included realistic effects such as 
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drag terms arising from the viscosity of the cytosohc fluid, and we have considered also only lateral 
(2-dimensional) diffusion in this work. 

It is seen that curvature effects modify considerably the solution of the diffusion problem. In 
the biological context, it is well known that depending upon the cell type, between 30-90 percent 
of all integral proteins in the plasma membrane are freely mobile and among these, the lateral 
diffusion rate of a protein in an intact membrane is around 10-30 times lower than that of the 
same protein embedded in synthetic liposomes [1,2,9]. It has been suggested that this could be 
because the mobility of the proteins might be hampered by interactions with the rigid submembrane 
cytoskeleton. 

In arriving at those diffusion rates, these authors considered only the normal form of Pick's law for 
flat space. Use of the correct form of Pick's law taking into account the varying curvatures of the 
membranes on which the protein molecules diffuse must be made when seeking to explain through 
such theories. 

Discussion 

We have shown how the curvature of the surface through which molecules diffuse modify the usual 
form of Pick's law and the relation between the concentration of the diffusing molecules, the diffusion 
constant and the time. 

Many intra-membrane particles are electrically charged and when they are subjected to an external 
electric field, move from their original random distribution to a more ordered distribution. It would 
be interesting to see the effect of external electromagnetic fields on molecules which are electrically 
charged, diffusing on curved surfaces. 

Acknowledgement 

I would like to acknowledge support from the Jawaharlal Nehru Centre for Advanced Scientific Re- 
search, Bangalore, during the course of this work. 



9 



References 

1. F.Zhang, G.M.Lee & K.Jacobson, Bioessays 15, 579 (1993) and references therein. 

2. K.Jacobson, A.Ishihara & R.Inman, Ann.Rev.Physiol.49, 163 (1987). 

3. A.E.Sowers & C.R.Hackenbrock, Biochim.Biophys.Acta 821, 85 (1985). 

4. P.Gilkey, J.Diff.Geom.lO, 601 (1975). 

5. B.S.DeWitt, Dynamical Theory of Groups and Fields, (Gordon & Breach, N.York, 1965). 

6. L.Parker and D.J.Toms, Phys.Rev.D31, 953 (1985). 
L.Parker and D.J.Toms, Phys.Rev.D31, 3424 (1985). 

7. L.Parker in Recent Developments in Gravitation - Cargese 1978 Lectures, eds. M.Levy and 
S.Deser (Plenum Press, New York, 1979). 

8. G.Gomppcr & J. Goes, in Annual Reviews of Computational Physics II -ed. D.Stauffer, 101- 
136. (World Scientific, Singapore (1995)). 

9. H.Lodish, D.Baltimore, A. Berk, S.Lawrence Zipursky, P.Matsudaira & J.Darnell, Molecular 
Cell Biology, (Scientific American Books, N.York, 1995) (3rd edition). 



10 



